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This paper describes the use of a Volume of Fluid (VOF) multiphase model to simulate 
the water injected from a rainbird nozzle used in the sound suppression system during 
launch. The simulations help determine the projectile motion for different water flow rates 
employed at the pad, as it is critical to know if water will splash on the first-stage rocket 
engine during liftoff. 


I. Introduction 

During the launch of the Space Shuttle several water suppression systems are used. These consist of the 
crest water nozzles mounted on the apex of the flame defector, Ignition Overpressure (IOP) water nozzles 
placed around the exhaust hole, and rainbird water nozzles installed above the launch deck (Figure 1). The 
rainbird water is specifically designed to suppress launch acoustics that can reflect from the flat surfaces of 
the deck. The Shuttle rainbirds were located on the North deck of the mobile launcher platform and were 
activated at T-0, so they didn’t really pose a threat to Shuttle. 

For the new Space Launch System (SLS) vehicle the operation is similar, regardless of the new mobile 
launcher and new engine configurations. The goal of the rainbird nozzle system remains the same, which 
is sound suppression (SS), and the rocket engines still cannot get wet. However, the rearrangement of the 
rainbird water system for the SLS mobile launcher, as shown in Figure 2, makes the rainbirds closer to the 
first-stage rocket engines, which are located above the exhaust hole and not shown in the Figure. The close 
proximity of the rainbird nozzle system could potentially cause vehicle wetting during liftoff. 

This study is to complement the one-dimensional projectile calculation 1 by considering the VOF multi- 
phase. The objective of this study is to determine if the Shuttle- heritage rainbird system can be used for 
the SLS without water splashing on its rocket engines. 

II. Modeling 

The problem under study involves free surface with sharp interface between two fluids. All modeling 
reported in this paper was performed using OpenFOAM 2.3.0, an open source CFD software package devel- 
oped by OpenCFD Ltd at ESI Group. 2 4 Available solvers for problems involving two incompressible fluids 
with an interface are either laminar, Reynolds- Average Navier-Stokes (RANS) turbulence, and Large Eddy 
Simulation (LES). Only laminar and RANS solutions are reported in this study. 

The problem considers the transient injection of water into an open space filled with air as shown in 
Figure 3. For simplicity we will consider only a two-dimensional case by assuming that the geometry has 
an infinite length in the z-direction. The rainbird nozzle with an exit diameter of 2.4 ft is canted at a 10 
degrees up and filled with water. At time zero, the water will be injected at an inlet velocity corresponding 
with a pre-determined volume flow rate. Table 1 lists the volume flow rates considered in the simulations 
for nomnial and abort operations. 

The throw distances based 1-D projectile for the North corner, North center, and South rainbird nozzles 
are shown in Figures 4, 5, and 6, respectively. The water volume flow for the nominal operation is based on 
a maximum rainbird flow during T-20 to T+20 seconds, and the abort case is the same as T-zero. 
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Rainbird Nozzle System 

Operation 

North East/ North West 

North Center 

South East/ South West 

Nominal 

72,491 

85,148 

115,065 

Abort 

35,783 

42,031 

56,798 


Table 1. Water Flow Rates in GPM 


III. Governing Equations 

In the VOF method the same continuity and momentum equations are solved for both gas and liquid 
phases. The physical properties of one fluid are calculated as weighted averages based on the volume fraction 
of the two fluids in one cell. The conservation form of continuity and momentum equations are given by: 

v-u = o 


dpU 

dt 


+ V • pUU — V • /iVU — pg = — Vp — S s 


where S s is the surface tension which takes place only at the free surfaces. 


The volume of fluid in a cell is computed as V = aV ce u , where V ce u is the volume of a computational cell 
and a is the fluid fraction in a cell. If the cell is completely filled with water then a = 1 and if it is filled 
with air then its value should be 0. At the interface the value of a is between 0 and 1. The phase fraction 
a can be computed from a separate transport equation that is given by: 


^ + V • aU = 0 
at 

A convective term is added to take into account of the compression of the surface, leading to a counter- 
gradient transport equation as follows: 


da 

~dt 


+ VmU + V-a(l- a)U c = 0 


where U c is a velocity field suitable to compress the interface. 


The density and viscosity in the domain are calculated as a weighted average of the volume fraction of 
the two fluids, cq as follows: 

p = api + (1 — a)p 2 
p = api + (1 - a)p>2 

The surface tension S s on the RHS of the momentum equation is given by: 


S a = aKVa 

K is the curvature of the interface and is calculated by: 


K = —V • h 

where n is the normal unit vector based on the volume fraction a: 

Va 
n = , 

|Va| 
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IV. Discussion of Results 


Figure 7 shows the structured mesh generated using blockMesh utility available in the OpenFOAM 
package. The same topology is used for both nominal and abort operations, except that the computational 
width is shorter in the abort case because of the shorter throw distance. The grading of the mesh can be 
easily modified to resolve the boundary layers close to the walls. However, this is not an important feature 
since the focus is on the actual distance traveled by the water jet. 

Figures 8-10 show the volume fraction of water (a) after 5 sec from the start of the injection for the 
higher volume flow rate associated with nominal operation; the red color corresponds to a=l, which is water 
and blue corresponding to a= 0, which is air. It can be seen that water stream has traveled 11m, 13m, and 
18.5m for the North corner, North center, and South rainbirds, respectively. These values closely match 
those depicted by the 1-D projectile calculation (11.9m, 14.6m, and 21.4m). The slight over predict ion by 
the projectile analysis is due to the fact that it neglected the air resistance effect, which is accounted for in 
the multi-phase modeling. 

Figures 11-13 show the velocity field for three rainbird nozzles during nominal operation. It is obvious 
that these simultions exhibit two different flow patterns. While the distance traveled can be approximated 
using these velocity contours and the results are in agreement with those obtained from the volume fraction 
contours, the high flow rate shows a less accumulation of water on the deck. For the lower volume flow rate, 
the recirculation zone is smaller and the flow is more diffusive at the impingement location. 

The contours of water volume fraction and velocity magnitude for the abort case can be found in Figures 
14-19. Similar conclusions can be drawn regarding the flow behaviors. Again, the throw distances are 
slightly overpredicted by the projectile calculation compared to multi-phase modeling; the distance traveled 
as predited by the VOF simulations are 5.2m, 6m, and 8m for North corner, North center, and South 
rainbirds, compared to 5.3m, 6.4m, and 8.9m from Figures 4-6. It should be noted that the laminar flow 
solution for the abort case is very dispersive. The turbulent diffusion solution was obtained using a k-e 
turbulence model and presented in the abort case. The throw distances for all rainbird nozzles operated 
under nominal and abort conditions are summarized in Table 2. 



Throw Distance (Nominal) 

Throw Distance (Abort) 

Analysis 

North East/ 

North 

South East / 

North East/ 

North 

South East / 

Model 

North West 

Center 

South West 

North West 

Center 

South West 

1-D Projectile 

11.9 m 

14.6 m 

21.4 m 

5.3 m 

6.4 m 

8.9 m 

Multi-phase 

11.0 m 

13.0 m 

18.5m 

5.2 m 

6.0 m 

8.0 m 


Table 2. Comparison of Throw Distances 

Both projectile and VOF predictions indicate that in the abort case the water jet stream from a rainbird 
nozzle will not reach the rocket engine. Likewise, during nominal operation it might appear that water 
jet from the rainbird operated at higher flow rates might reach the core stage. However, the vehicle starts 
to move before the water jet reaches fully-developed flow; therefore, water impingement on the vehicle is 
unlikely. A detailed three-dimensional VOF or smoothed particle hydrodynamics (SPH) analysis 5-8 could 
be helpful in studying such transient phenomenon. 

V. Conclusion 

The 2-D VOF simulation has been used to determine the throw distances by the rainbird water nozzle 
system. It provides a less conservative prediction as compared to 1-D projectile analysis, thus allowing 
design margins to be reduced. In general, there is good agreement between the two analyses for the short 
throw distances. The discrepancies are greater for long throw distances, in which air resistance and drag 
dominated flows are more significant. The IOP/SS engineers have several design options in case the rainbird 
water injection system could wet rocket engines during launch. These options include lowering the rainbird 
nozzle, adjusting the flow rates, or changing its injection angle. Further study using three-dimensional VOF 
or SPH might be required to determine if Shuttle-derived rainbirds can meet the requirement of sound 
suppression while keeping the SLS rocket engines and other sensitive ground support elements dry during 
nominal launch and pad abort operations. 
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Figure 1. Shuttle Water Sound Suppression System. 
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Figure 2. SLS Rainbird Nozzle System. 
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Figure 3. Geometry Setup. 
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Figure 4. 1-D Projectile Prediction (North Corner). 
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Figure 5. 1-D Projectile Prediction (North Centerr). 
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Figure 6. 1-D Projectile Prediction (South). 
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Figure 7. Computatioal Mesh. 
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Figure 8. Water Volume Fraction for North Corner Rainbird (Nominal - 72,491 GPM). 
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Figure 9. Water Volume Fraction for North Center Rainbird (Nominal - 85,148 GPM). 
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Figure 10. Water Volume Fraction for South Rainbird (Nominal - 115,065 GPM). 
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Figure 11. Velocity Magnitude for North Corner Rainbird (Nominal - 72,491 GPM). 
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Figure 12. Velocity Magnitude for North Center Rainbird (Nominal - 85,148 GPM). 
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Figure 13. Velocity Magnitude for South Rainbird (Nominal - 115,065 GPM). 
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Figure 14. Water Volume Fraction for North Corner Rainbird (Abort - 35,783 GPM). 
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Figure 15. Water Volume Fraction for North Center Rainbird (Abort - 42,031 GPM). 
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Figure 16. Water Volume Fraction for South Rainbird (Abort - 56,798 GPM). 
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Figure 17. Velocity Magnitude for North Corner Rainbird (Abort -35,783 GPM). 
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Figure 18. Velocity Magnitude for North Center Rainbird (Abort - 42,031 GPM). 
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Figure 19. 
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Velocity Magnitude for South Rainbird (Abort - 56,798 GPM). 
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